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Abstract 

It has recently been observed that the SIR distributions of a variety of cellular network models and transmission 
techniques look very similar in shape. As a result, they are well approximated by a simple horizontal shift (or 
gain) of the distribution of the most tractable model, the Poisson point process (PPP). To study and explain this 
behavior, this paper focuses on general single-tier network models with nearest-base station association and studies 
the asymptotic gain both at 0 and at infinity. 

We show that the gain at 0 is determined by the so-called mean interference-to-signal ratio (MISR) between 
the PPP and the network model under consideration, while the gain at infinity is determined by the expected 
fading-to-interference ratio (EFIR). 

The analysis of the MISR is based on a novel type of point process, the so-called relative distance process, 
which is a one-dimensional point process on the unit interval [0,1] that fully determines the SIR. A comparison of 
the gains at 0 and infinity shows that the gain at 0 indeed provides an excellent approximation for the entire SIR 
distribution. Moreover, the gain is mostly a function of the network geometry and barely depends on the path loss 
exponent and the fading. The results are illustrated using several examples of repulsive point processes. 

Index Terms 

Cellular networks, stochastic geometry, signal-to-interference ratio, Poisson point processes. 


A. Motivation 


I. Introduction 


The distribution of the signal-to-interference ratio (SIR) is a key quantity in the analysis and design of interference- 
limited wireless systems. Here we focus on general single-tier cellular networks where users are connected to the 
strongest (nearest) base station (BS). Let <I> C be a point process representing the locations of the BSs and let 
xq C ‘J’ be the serving BS of the typical user at the origin, i.e., define xq = arg min{x G ||x||}. Assuming all 
BSs transmit at the same power level, the downlink SIR is given by 


5IP A ^ f^XQ^(xo) 

I ExG^\{xo}hxi{x)' 


( 1 ) 


where (hx) are iid random variables representing the fading and i is the path loss law. The complementary cumulative 
distribution (ccdf) of the SIR is 


FsirW=F(SIR>0). 


( 2 ) 


Under the SIR threshold model for reception, the ccdf of the SIR can also be interpreted as the success probability 
of a transmission, i.e., Ps{0) = -^sir(^)- 

In the case where <1> is a homogeneous Poisson point process (PPP), Rayleigh fading, and i{x) = ||x||““, the 
success probability was determined in ||2l. It can be expressed in terms of the Gaussian hypergeometric function 
22^1 as El 


Ps,ppp(0) = 


1 


2Fi(1,-5;1-5;-0)’ 


(3) 
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where 5 = 2!a. For a = 4, remarkably, this simplifies to 

Ps,ppp(^') = -7=-7=- 

1 + y/O arctan y/O 

In im, it is shown that the same expression holds for the homogeneous independent Poisson (HIP) model, where 
the different tiers in a heterogeneous cellular network form independent homogeneous PPPs. For all other cases, 
the success probability is intractable or can at best be expressed using combinations of infinite sums and integrals. 
Hence there is a critical need for techniques that yield good approximations of the SIR distribution for non-Poisson 
networks. 


B. Asymptotic SIR gains and the MISR 

It has recently been observed in Q, Q that the SIR ccdfs for different point processes and transmission techniques 
(e.g., BS cooperation or silencing) appear to be merely horizontally shifted versions of each other (in dB), as long 
as their diversity gain is the same. 

Consequently, the success probability of a network model can be accurately approximated by that of a reference 
network model by scaling the threshold 6 by this SIR gain factor (or shift in dB) G, i.e., 


Ps{0) ~Ps,ref(6'/G)- 


Formally, the horizontal gap at target probability p is defined as 

pe (0,1), 


Op(p) ^ 




SIRref 


(P)’ 


(4) 


where Fg|p is fhe inverse of fhe ccdf of fhe SIR and p is the success probability where the gap is measured. It is 
often convenient to consider the gap as a function of 6, defined as 

G{e) 4 Gp(p,,ref(0)) = ^5IRfe.ref(g)) ^ 

Due to its tractability, the PPP is a sensible choice as the reference modeiH 

If the shift G{6) is indeed approximately a constant, i.e., G{9) G, then G can be determined by evaluating 
G{9) for an arbitrary value of 9. As shown in ||6l, the limit of G{9) as 0 —)• 0 is relatively easy to calculate. 
Here we focus in addition on the positive limit 9 ^ oo and compare the two asymptotic gains to demonstrate the 
effectiveness of the idea of horizontally shifting SIR distributions by a constant. 

So the main focus of this paper are the asymptotic gains relative to the PPP, defined as follows. 


Definition 1 (Asympfolic gains relative fo PPP). The asymptotic gains (whenever the limits exist) Gq and Goo are 
defined as 

Go = lim G{9)] Goo = lim G{9), (6) 

0 —^0 0 —^oo 

where the PPP is used as the reference model. 


C. Prior work 

Some insights on Gq are available from prior work. In @ it is shown that for Rayleigh fading, Gq is closely 
connected to the mean interference-to-signal ratio (MISR). The MISR is the mean of the interference-to-(average)- 
signal ratio ISR, defined as 

ISR = 


I 


Eh(5) ^ 


'This is why the method of approximating an SIR distribution by a shifted version of the PPP’s SIR distribution is called ASAPPP— 
“Approximate SIR analysis based on the PPP” (7). 
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Fig. 1. The SIR distributions for the PPP (solid) and the triangular lattice (dashed) for a = 4 and the lower bound (which is asymptotically 
tight) for the PPP (dash-dotted). The horizontal gap between the SIR distributions of the PPP and the triangular lattice is 3.4 dB for 
a wide range of 9 values. The shaded band indicates the region in which the SIR distributions for all stationary point process fall that are 
more regular than the PPP. 


where E/j(5) = E(S | rh) is the mean received signal power averaged only over the fading. Not unexpectedly, the 
calculation of the MISR for the PPP is relatively straightforward and yields MISRppp = 2/(a — 2) 10 Eqn. (8)J1. 
Since P(SIR > 9) = P{h > 6 ISR) = E(P(/i > 6 ISR | ISR)), the success prohahility can in general he expressed as 

Ps{9)=EFh{e\SR), (7) 

where Fh is the ccdf of the fading random variables. For Rayleigh fading, Fh{x) = and thus Ps{9) r\j 1-0 MISR, 
0 —)• 0, resulting in 

MISRppp 2 1 

° ^ MISR ^ a-2 MISR 

and 

Ps(0) ~ Ps,ppp(^'/G'o), 0 —)• 0. 

So, asymptotically, shifting the ccdf of the SIR distrihution of the PPP is exact. 

An example is shown in Fig. [T] where 0 = 4, which results in MISRppp = 1, while for the triangular lattice 
MISRtri = 0.457. Hence the horizontal shift is MISRppp/MISRtri = 3.4 dB. For Rayleigh fading, we also have 
the relationship Ps{0) = 71|5r( 0) > g-^'MiSR Jensen’s inequality, also shown in the figure. Here is a lower 
hound with asymptotic equality. 

In lUl, the authors considered coherent and non-coherent joint transmission for the HIP model and derived 
expressions for the SIR distrihution. The diversity gain and the asymptotic pre-constants as 0 —)■ 0 are also derived. 
In Q, the henefits of BS silencing (inter-cell interference coordination) and re-transmissions (intra-cell diversity) 
in Poisson networks with Rayleigh fading are studied. For 0 —)• 0, it is shown that Ps{9) ~ 1 — akO when the k — 1 
strongest interfering BSs are silenced, while Ps{9) ~ 1 — for intra-cell diversity with m transmissions. For 

0 -)■ CX), Ps{0) AkO ^ and ps(0) r\j BmO ^ for BS silencing and retransmissions, respectively. The constants a^, 
bm, Ak, and Bm are also determined. Fastly, fS) Thm. 2] gives an expression for the limit lime^oo 9^ps{6) for the 
PPP and the Ginihre point process (GPP) with Rayleigh fading. For the GPP, it consists of a double integral with 
an infinite product. 

In 121, the authors consider a Poisson model for the BSs and define a new point process termed signal-to-total- 
interference-and-noise ratio (STINR) process^ They obtain the moment measures of the new process and use them 


different derivation of this result will be given in Thm. [T] in Sec. IV in this paper. 

^What is meant by “total interference” is actually the total received power (including the desired signal power). 
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to express the probability that the user is covered by k BSs. In our work, we consider a different map of the original 
point process based on relative distances, which results in simplified moment measures for the PPP and permits 
generalizations to other point process models for the base stations. 


D. Contributions 


This paper makes the following contributions: 

• We define the relative distance process (RDP), which is the relevant point process for cellular networks with 
nearest-BS association, and derive some of its pertinent properties, in particular the probability generating 
functional (PGFL). 

• We introduce the generalized MISR, defined as MISR„ = (E(ISR"^))^/”, which is applicable fo general fading 
models, and give an explicif expression and fighf bounds for fhe PPP 

• We provide some evidence why fhe gain Go is insensitive to the path loss exponent a and the fading statistics. 

• We show that for all stationary point process models and any type of fading, the tail of the SIR distribution 
always scales as 6~^, i.e., we have Ps{0) ~ c6~^, 6 —)■ oo, where the constant c captures the effects of the 
network geometry and fading. The asymptotic gain follows as 


and we have 



Ps{0) ^ Ps,ppp{6/Goo), 9^00. 


( 8 ) 


• We introduce the expected fading-to-interference ratio (EFIR) and show that the constant c is related to the 
EFIR by c = EFIR'^. Consequently, Goo is given by the ratio of the EFIR of the general point process under 
consideration and the EFIR of the PPP. 


II. System Model 

The base station locations are modeled as a simple stationary point process d* C M^. Without loss of generality, 
we assume that the typical user is located at the origin o. The path loss between the typical user and a BS at x E 
is given by i{x) = ||x||““, a > 2. Let denote the ccdf of the iid fading random variables, which are assumed 
to have mean 1. 

We assume nearest-BS association, wherein a user is served by the closest BS. Let xq denote the closest BS to 
the typical user at the origin and define R = ||xo|| and <!>■ = <I> \ {xq}. Wifh fhe nearest-BS association rule, the 
downlink SIR ([T]) of the typical user can be expressed as 


SIR 


hR-^ 


Further notation: b{o, r) denotes the open disk of radius r at o, and b{o, rY 


(9) 

= \ 6(o, r) is its complement. 


III. The Relative Distance Process 

In this section, we introduce a new point process that is a transformation of the original point process and 
helps in the analysis of the interference-to-signal ratio. 


A. Definition 

Prom ([T]), the MISR is defined as 

MISR = E f • (10) 

V ^(^o) / V J 

The firsf expecfafion is faken over <1> and h, while fhe second one is only over d> since E(/i) = 1. Since £(x) only 
depends on ||x||, it is apparent that the MISR is determined by the relative distances of the interfering and serving 
BSs. Accordingly, we introduce a new point process on the unit interval (0,1) that captures only these relative 
distances. 
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Definition 2 (Relative distance process (RDP)). For a simple stationary point process let xq = arg min{x G 
<P: ||x||}. The relative distance process (RDP) is defined as 

7^4{xGd>\{xo}: ||xo||/||x||}c(0,l). 

Using the RDP, the ISR can he expressed as 


and, since IE((iy) = 1, the MISR is 


ISR = hyV^, 



For the stationary PPP, the cdfs of the elements vj. = ||xo||/||xfc|| of TZ are F,y^{x) = 1 — (1 — x"^)^, x G [0,1], 
as given in |j6|. Summing the densities dF),j,(x)/dx over A: G N yields the intensity measure A(dr) = 2r~^dr. It 
follows that the mean measure A([r, 1)) = E7^([r, 1)) = — 1, 0 < r < 1. The fact that the mean measure 

diverges near 0 is consistent with the fact that TZ is not locally finite on intervals [0, e). 

Generally, since the MISR only depends on the relative distances, the gain Go does not depend on the base 
station density. 


B. RDP of the PPP 

The success prohahility for Rayleigh fading is given hy the Laplace transform of the ISR: 

Ps( 0 ) = Ee-"'^^ = E =E J] —( 11 ) 

yen yen ^ ^ 

This RDP-hased formulation has the advantage that it circumvents the usual two-step procedure, where first the 
conditional success prohahility given the distance to the serving base station R is calculated and then an expectation 
with respect to R is taken. 

Since (ITTI) has the form of a PGFL, we first calculate the PGFL of the RDP generated by a PPP. 


Lemma 1. When d? is a PPP, the probability generating functional of the RDP is given by 

Gnif] = E n fix) = -^- - -, 

Ln l + 2/o\l-/(x))x-3dx’ 

for functions f: [0,1] i-G [0,1] such that the integral in the denominator of (1121) is finite. 
Proof: The PGFL G-jilf] can be calculated as 


Gn[f] = E 


= E 


n/(^ 

.xen 

n / 


|xo| 


X 


E/jexp 27rA J a — / ^ — 1 1 do 

OO / OO 


R 


(b) 


= A J 27rrexp ^ —27rA J o (^1 — / |^- ) ) do | e dr 


CXJ 


= Ay 27rrexp 27rr^A J y ) ) dy | e dr 

1 


1 + 2/ry(l - /(l/ 2 /))dy’ 


( 12 ) 
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where (a) is obtained from the PGFL of the PPP, which, for a general PPP 'P with intensity measure A, is given 
by 



(13) 


Writing the PGFL in polar form and conditioning on the distance to the nearest neighbor R= ||xo|| yields (o). In 
{h), we de-condition on R using the the nearest-neighbor distribution of the PPP Using the substitution y~^ —)■ x, 


we obtain the final result. 


When f{x) = 1/(1 -I- 6x°‘) (see (fTTIl l. we retrieve the result in (l3]l for Poisson cellular networks with Rayleigh 
fading. 

It may be suspected that the RDP of a PPP is itself a (non-stationary) PPP on [0,1]. It is easily seen that this 
is not the case. Let 'P be a PPP on [0,1] with the same intensity function as TZ, i.e., A(dr) = 2r~^dr. If TZ was 
a PPP, the success probability for Rayleigh fading would follow from the PGFL of 'P (specializing ([T3]) to a PPP 
on [0,1]) as 



[/] = exp 


and, in turn, the success probability would be given by 



(14) 


instead of (l3]l. 


However, assuming TZ to be Poisson yields an approximation of the success probability, with asymptotic equality 
as 0 —)■ 0. The “Poisson approximation” (fT4l) is related to the actual value © as 



This holds due to the identity 


69 2 F 1 {1,1-d-,2-6;-9) 
1-6 


= 2Fi(1,-5;1-,5;-0)-1. 


Rewriting and expanding, we have 



1 





Hence, only considering the dominant first term in the denominator as 0 —)■ 0, we obtain Ps{9) rsj Ps{9), 9^0. 

The fact that Ps{9) < Ps{9) for 0 > 0 is an indication that the higher moment densities of the RDP are larger 
than those of the PPP. This is indeed the case, as the following calculation of the moment densities shows. 

Lemma 2. When <P is a PPP, the moment densities of the RDP are given by 


n 




( 15 ) 


Proof: First we obtain the factorial moment measures. We use the simplified notatior0 

a^'^\ti,t2, . . . ,fn) = 1) X (f2, 1) X ... X {tn, 1)), 0 < < 1. 


The factorial moment measures are defined as 





''Here the intervals are chosen as {t, 1) for t > 0 since the RDP is not locally finite on [0, e]. 
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where indicates that the sum is taken over n-tuples of distinct points. The moment measures are related to the 
PGFL as Uni p. 116] 

• • • ,fn) = ~ “ ®2lp2,l) - • • • - (16) 

evaluated at si = S 2 = • • • = Sn = 0. Using Lemma [T] we obtain 

Gtz[1 - - S2l(t„l) - ... - Snl(t„,l)] = ^ ' 

Differentiating with respect to Si and setting si = S 2 = • • • = Sn = 0, we have 

a^’^)(fl,f2, ■ ■ ■ Un) =^^11 

j=i V i / 


The moment densities follow from differentiation, noting that fj denotes the start of the interval, which causes a 
sign change since increasing fj decreases the measure. ■ 

So the product densities are a factor n! larger than they would he if TZ was a PPP. This implies, interestingly, 
that the pair correlation function ifTTl Def. 6.6] of the RDP of the PPP is g{x,y) = 2, \/x,y G (0,1). 

The moment densities of the RDP provide an alternative way to obtain the success probability for the PPP: 

1 




yen 


yen 


1 + 


= l-'^u{9,y) + —E ^ u{9,yi)u{e,y2) + ... +^ Y[u{9,yi) + ..., 

yen yi,y2en yi,...,yneni=i 

where i'{9,y) = . From the definition of the moment densities, we have 


n=0 

oo 

E 

n=0 


7^ 


yi,...,ynen *=1 

(- 1 )- 


n! 


'[ 0 , 1 ]" 




Ki=l 


Using Lemma O we have 


OO / \ ^ 

Ps(0) = , 

^0 vm J 


1 - 5 

^9^{-ir - 6--9) - 1) 

1 


n=0 

oo 


n=0 


2Fi(1,-5;1-5;-0)’ 

which equals the success probability given in Q- 


(18) 


C. RDP of a stationary point process 

We now characterize the PGFL of the RDP generated by a stationary point process. Let f{R, <!>’) be a positive 
function of the distance R = ||xo|| and the point process d>' = <1> \ {xq}- The average E[f{R, <!>■)] can in principle 
be evaluated using the joint distribution of R and which is, however, known only for a few point processes. 
Thus we introduce an alternative representation of f{R, <h') that is easier to work with. 











The indicator variable l(<h(6(o, ||x||)) = 0), x G equals one only when x = xq. Hence it follows that 

f{R, CD') = /(||x||, CD \ {x})l(cD(6(o, ||x||)) = 0). (19) 




This representation of f{R, cD') permits the computation of the expectation of f{R, cD') using the Campbell-Mecke 
theorem ifTTl Thm. 8.2]. We use the above idea in the next lemma to obtain the PGFL of a general RDP. 


Lemma 3. The PGFL of the RDP generated by a stationary point process cD is given by 


Gn[f] = \ Q[ 


f 


■ +x\\J 




(• + X G b{o, ||x||)^) 


dx, 


( 20 ) 


where is the PGFL of the point process cD with respect to the reduced Palm measure. 
Proof: Using the representation in ([19]), we obtain 


Gn[f]=^Y. n f 

A / KUf 

/ 


(b 


(a) 


= A 




\\y + x\ 
||x|| 

l|y + a:|| 


l(T>(6(o, ||x||)) =0) 

1(^>(6(—X, ||x||)) = 0)dx 
l{y G b{-x, ||x||)'')dx 


/ 


= X Qr 


f 


X 


+x|| 


!(• + X G b{o, ||x||)^) 


dx, 


where (a) follows from the Campbell-Mecke theorem. ■ 

If cD is also rotationally invariant {i.e., motion-invariant), the reduced Palm measure is also rotationally invariant 
and hence 

CXD 

^ ^ r 


Gn[f] = A 27 r I rG[ 
0 


/ 


• +^'11 


!(• -I- r G b{o,ry) 


dr, 


( 21 ) 


where r in any vector addition should be interpreted as (r, 0). 
Remark. Taking cD to be a PPP, (|2TI) reduces to ([12]) since 


Gnlf] = 2vrA [ rexp [ - A / 
Jo \ 

OC 

■I 


1 -/ 


= 27rA / r exp ( — Avrr^ — A 


b{o,r)^ 


|x-h (r,0)|| 


1-/ ^ 
X 


l(x -I- (r, 0) G 6(0, rY) 


dxj dr 


dx dr 


= 27rA J r exp A7rr^ — 27rA J y[l — f{r/y)]dyj dr. 

The last expression equals the second-to-last line in the proof of Lemma [T] 

Similar to the case of PPP, we now obtain the moment measures of the RDP of a general stationary point process. 

Lemma 4. The factorial moment measures of the RDP TZ generated by a stationary point process T* are 

P n 

= J Y{f{yo,yi,ti)^y„^y,^,„^yJg{^,yo)]p^^^^\yo,yi,---,yn)dyo---dyn, (22) 

Rn+l *=1 

where p('D,x) = flyg'!' l(lly|l ^ ll^ll) f{x,y,t) = Ipq) (||f{[)- product densities are 

■ ■ ■ fn) =n!2” 1 Yn{tl,... Rn), (23) 


•^ 2=1 


















where 


' / 

n\ 2 ^ J 

II2/0II 

1 

R2 


[0,27r] 


fl/O, 

(hoW 


1 G ’ 


l), ..., d^Ji.. .dvJndyo- (24) 


Proof: As before, we use the relationship (fT^ . While the result can be obtained from the PGFL in Lemma [H 
it is easier to begin with the definition of the PGFL. We have 

0r|/i = e 5; n 

We are interested in the derivative of the PGFL with the function f{z) = 1 — Y^^=i So we have 

d d 




dsi ds, 


-E^g{<^>\{x},x) JJ f 1 - ^Sil{t.,i) 

'■ IG# 3/G‘I’\{2^} ' *=t 


||X|| 

M 


evaluated at si = S 2 = ... = Sn = 0. Expanding the inner product over the summation we obtain an infinite 
polynomial in the powers of si,... ,Sn and their products. We observe that the only term that contributes to the 
derivative in a non-zero manner is the siS 2 • • • Sn term. Its coefficient equals 


T = (-l)"Er^5r(4> \ {x},x) ^ f{x,yi,ti)...f{x,yn,tn)\ 

xG$ yi,...y„G^\{x} 

Combining the summations, 

r = (-i)’^Er ^ f{yo,yi,ti)- ■■■■ f{yo,yn,tn)g{^\{yo},yo)\ 

'' ->/_ n, zrrIS ' 


' yo,yi,---yn&^ 

Since f{yo,yi,ti) / 1 implies ||7/i|| > ||yo||, 

7^ 


r=(-l)”Ef f{yo,yi,ti)-...■ f{yo,yn,tn)g{^\{yo,yi,...,yn},yo)\ 

^ V0,Vl,---VnG$ '' 


and, using ifT^ Thm. 1], 


T = {-1)"- I Y{f{yo,yi,ti)E!y^ y^ yJg{^,yo)] ■ pY^\yo,yi, ■ ■ ■ ,yn)dyo ■ ■ ■ dyn, 


(25) 


2 = 1 


and the result (l22l) follows. For the product densities, we convert the variables Xi into polar coordinates {ri,cpi), 
which yields 

llVoli/ti llvoll/L 

llvoll llvoll [ 0 , 2 ^" 

• {rn,‘Pn))d(pi .. .dy^ndri... dr^dyo- 

Then differentiating using the Leibniz rule with respect to fi,..., we obtain 


p5^^(tl,t2,...,fn) = (nti /lli/of" f 
\i=i J 


[0,27r]’^ 


A.(—. 


(n-l-1) 

P\, 2/0, 


Wl 


, <^1 / , • • •, 


II2/0II 


ifn difi. ..dipndyQ, 
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which equals (|2^ . ■ 

As in (fT^ . the moment densities of the RDP generated hy a stationary point process can he used to compute its 
corresponding success prohahility. 

IV. The Mean Interference-to-Signal Ratio (MISR) and the Gain at 0 

In this section, we introduce and analyze the MISR, including its generalized version, and apply it to derive a 
simple asymptotic expression of the SIR distribution near 0 using the gain Gq. We also give some insight why Go 
harely depends on the path loss exponent a and the fading statistics. 


A. The MISR for general point processes 

The first result gives an expression for the MISR for a general point process. 
Theorem 1. The MISR of a motion-invariant point process <I> is given by 

MISR = 2/ 

Jo 

where /3i (t) is given in (1241) in Lemma H] 


Proof: Using the RDP, the MISR can he expressed as 

MISR = E^y“ 

yen 

= I 

Jo 

= 2 [ 

Jo 

where (a) follows from Lemma |4l ■ 

When is a PPP, from Slivnyak’s theorem and the fact that p^^'i = A^, we have = 1 and hence MISR = 
2/(a-2). 


B. The Generalized MISR 

Definition 3 (Generalized MISR). The generalized MISR with parameter n is defined as 

MISRn = (E(ISR’"))^/”. 

If there is a danger of confusion, we call MISR = MISRi the standard MISR. 

The generalized MISR can he obtained by taking the corresponding derivative of the Laplace transform E(e“^ 
at s = 0. In case of the PPP with Rayleigh fading, the Laplace transform is known and equals the success probability 
dD, thus 

misr;^ = E(ISR") = (-1)-(26) 

For general fading, the Laplace transform is not known, but we can still calculate the derivative at s = 0, as the 
following result for the PPP with general fading shows. 


Theorem 2 (Generalized MISR and lower bound for PPP). For a Poisson cellular network with arbitrary fading, 

5 6E{h^-^+^) 


E(ISR") = 


fc=i 


l-(5’ 


n — k 1 — 6 


where Bn,k ^re the (incomplete) Bell polynomials. For n > 1, the generalized MISR is lower bounded as 


MISR„ > 


1-5 


n! + 


5E{h^) 

n — 6 


l/n 


(27) 


(28) 
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For n = 2, equality holds, and for <5 —)• 0 and 5 —>• 1, the lower bound is asymptotically tight. 
Proof: We begin with the the Laplace transform of the ISR, given by 

yeti 

(a) 1 


1 + 2/-, 

where (a) follows from Lemma[T] Let f{s) = 1/(1+ s) and g{s) = 2 / — , (^1 — C 


-slSR^ 


d,. Then E(e 

f{g{s)). We are interested in the m-derivative of E(e“^'^^) with respect to s at s = 0, which can be computed 
using Faa di Bruno’s formula lIT^ as 

1 n 

-Jig{s))U=o = ^f^^H9m)BnA9'i0),g"i0),...,g^^-’^+^\0)), 


k=l 


where the (incomplete) Bell polynomials. We have ,(0) = 0, 

+ kl 


and 


/ oo 

which, when evaluated at s = 0, equals 

9^^>{0) = 

Combining everything, we have 


2(-l)^+^E(fe^) 
ka — 2 


- 

E(I5R") = (-1)"— 

n 




k=l 


2 2(-l)’^-^E(/i”-^+t) 

a — 2 ’ ’ {n — k + l)a — 2 


(29) 


From the definition of Bell polynomials it follows that all the terms are positive, hence the result (1271) follows from 
S = 2/q!. The lower bound is obtained by only considering the terms k = 1 and k = m in the sum (l29l ). The 
bound becomes tight as 5 —)■ 0 and as <5 —)■ 1 since the term k = 1 dominates the sum (l29l ) as J —>• 1 since it is the 
only term with a denominator (1 — 5)”, while the term k = n dominates as 5 —0 since it is the only one with a 
numerator 0(5). ■ 

Hence we have two simpler asymptotically tight bounds for the generalized MISR: 


MISRn > ( -E(/i" 
n 


l/n 


0 


MISR„> °V-\ = MISRi(n!)V", 5^1. 
1 — 0 


(30) 

(31) 


For Rayleigh fading, (l30l) yields MISR„ ~ (5r(n))^/’^, 5^0. 

Fig. |2] shows MISR„ for Rayleigh fading as a function of the path loss exponent. As can be observed, the term 
MISRi(n!)^/” is dominant for a < 4 even if the fading is severe (Rayleigh fading). For less severe fading, the term 
with E(h^) is less relevant; it only becomes dominant for unrealistically high path loss exponents (5 <C 1). 

The second moment of the ISR follows from (|27]) as 

E(ISR^) = 2 MISR? + 

2 — 0 
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Fig. 2. MISR„ for n € {1, 2, 5} for the PPP with Rayleigh fading as a function of the path loss exponent a. “Lower bound 1” is the 
simple bound in 113, which holds irrespective of the fading and is asymptotically tight as a J, 2, and “Lower bound 2” is the bound in l l28t . 
which is valid for n > 2 and is exact for n = 2. For a < 4, the two bounds are essentially identical. 


and the third moment is 


E(ISR^) = 6MISR? + 


6<52E(/i2) 

(1-5)(2-<5) 


+ 


(5E(/i3) 

3-<5 ■ 


Remarks. 

• Setting = 1 for all k retrieves the result in l[3l Prop. 3] on the pre-constant for m transmissions in a 

Poisson networks over Rayleigh fading. 

• An alternative way to derive the lower hound is as follows. Letting Uy = hyy‘^ for y G 7^, we expand ISR"^ as 


ISR'" = 


Un 





UxUz + 


'Xi ^Xn ? 


where the expression contains k sums. Ignoring all hut the first and last terms of the expansion, we obtain 

E(ISR") > E(/r"')— [ {xi---Xn)°‘p^^\xi,...,Xn)dxi---dxn 

m-6 7(0,i]n 

= E(/i’")^—+ n!MISR^ 
n — 0 

which equals the result in (l28l ). 

For Nakagami-m fading, MISR„ is decreasing with increasing m since the moments K{h^) are decreasing with 
m. As the lower hound MISRi(n!)^/’^ does not depend on the fading, MISR„ approaches a non-trivial limit as 

m ^ oo. 

Fig. [3] shows MISR„ as a function of n. The increase is almost linear in n. Indeed, as n —>• oo, MISR„ is 
proportional to n for the usually encountered path loss exponents, as the following corollary estahlishes. 
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Fig. 3. Generalized MISR per ^Tl\ for the PPP for Nakagami-m fading with m £ {1, 2,10}. 


Corollary 3. For the PPP with Rayleigh fading and a < 4, 

Tl Tl S 

MISRn ~ — MISRi =- r, n oo. 

e e 1 — 0 

Proof: For the PPP with Rayleigh fading and 5 > 1/2, it follows from (l27l) that 


MISR„ ~ 



n —>■ oo, 


since the dominant term in (|27]) for large n is the one with 5”/(l — 5)'^, which increases geometrically (or stays 
constant) with n for <5 > 1/2. For the factorial term, log((n!)^/”) ~ logn — 1, hence we obtain MISRn ~ 
giogn-i MISRi. ■ 

Remark. Using Stirling’s formula n! ~ y/2TTn{n/e)'^, this asymptotic result can he sharpened slightly. 


C. The gain Gq for general fading 

Equipped with the results from Theorem^ we can now discuss the gain Gq for general fading If Fh{x) ~ Cmx"^, 
X —)• 0, then, for 0 —)• 0, we have PsiO) ~ 1 — Cm ISR)™], hence 

(„,) _ /E(ISR^\ _ MISRn,,PPP 

° V E(ISR™) ) MISRn, ■ ^ 

The ASAPPP approximation follows as 

!>.(«) “pSppIS/gM. 

where P^^pp is the success prohahility for the PPP with fading parameter m, which is not known in closed-form. 
In ifTdl . the SIR ccdf for a Poisson cellular network when h is gamma distributed is discussed. However, we have 
the exact MISR^ from (ITtI) and the lower bound MISRn, ^ MISRi(m!)^/”^. 

For Nakagami-m fading, the pre-constant is Cn, = m™“^/r(m), and we have 


pS?ppW~l-CmE[(0ISR)”^] 

m-1 

<l--p^MISRim!0”^ 

1 [m) 

= 1 - MISRi(m0)”^, 


^By “general fading”, here we refer to a fading distribution that satisfies Fh{x) = ©(a:"*), a; —>■ 0, for arbitrary m £ N. 
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Fig. 4. Gain of selection combining of m transmissions in Rayleigh fading channels over a single transmission in Nakagami-m fading 
channels, for Poisson networks with a = 4. 


where indicates an upper bound with asymptotic equality. Adding the second term in the lower bound and 
noting that 




yields the slightly sharper result 


mS 

1-6 


+ 


r(2m) 


m — 6 T{rri)rrf‘ 


The gain for general fading is applicable to arbitrary transmission techniques that provide the same amount 
of diversity, not just to compare different base station deployments. As an example, we determine the gain from 
selection combining of the signals from m transmissions over Rayleigh fading channels with a single transmission 
over Nakagami-m fading channels, both for Poisson distributed base stations. The MISR for the selection combining 
scheme follows from 13] Prop. 3]. Fig. |4] shows that there is a very small gain from selection combining. 

Simulation results indicate that at least for moderate m, the scaling MISRm ~ MISRi(m!)^/™ holds for arbitrary 
motion-invariant point processes. This implies that which indicates that Go is insensitive to the fading 

statistics for small to moderate m. Next we show that the gain is also insensitive to the path loss exponent a. 


D. Insensitivity of the MISR to a 

Fig. [5] illustrates the densities of the square and triangular lattices relative to the PPP’s, which is Appp(r) = 2r~^, 
r G (0,1], as derived after Def. |2] Since the relative densities are roughly constant over the [0,1] interval, the gains 
do not depend strongly on a. Indeed, if the density of the RDP of a general point process could be expressed as 
A(r) = cAppp(r), we would have Go = 1/c irrespective of a. 

Another way to show the insensitivity of the gain to a is by exploring the asymptotic behavior of the MISR for 
general point processes given in Theorem [T] in the high-a regime. The result is the content of the next lemma. 

Lemma 5. For a motion-invariant point process <F, 

MISR((i) ~ (5/3i(l), (5^0, (33) 

/3i(l) = \ (yo,(||yo||,<Pi))d(/?idyo. 

[0,27r] 


where 
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Relative density of RDP for square lattice Relative density of RDP for triangular lattice 



(a) Relative intensity Asq(r)/Appp(r) of square lattice. (b) Relative intensity Atri(f)/Appp(r) of triangular lattice. 


Fig. 5. Relative intensity of square and triangular lattices. The straight line is the average relative intensity and corresponds to l/Go.sq 
and l/Go,tri, respectively. 


Proof: The MISR for a general point process is given by Theorem [T] as 


Jo Jo 


MISR = 2 

to Jo 

Using the Laplace asymptotic technique ifTSl Eq. 6.419], 

MISR(a) ~ 2a“^/3i(l), a ^ oo. 


(34) 


This shows that MISR for arbitrary point processes decays as 1/a, which implies Gq approaches a constant for 
large a (see Fig. Oa)). 


V. The Expected Fading-to-Interference ratio (EFIR) and the Gain at oo 

In this section, we define the expected fading-to-interference ratio (EFIR) and explore its connection to the gain 
Goo in ([D. We shall see that the EFIR plays a similar role for 6* —)• oo as the MISR does for 0 —)• 0. 


A. Definition and EFIR for PPP 


Definition 4 (Expected fading-to-interference ratio (EFIR)). For a point process <1>, let loo = I 

let h be a fading random variable independent of all (hx)- The expected fading-to-interference ratio 
defined as 


EFIR = 




x\\ " and 
(EFIR) is 


(35) 


where is the expectation with respect to the reduced Palm measure of ‘h. 


Here we use loo for the interference term, since the interference here is the total received power from all points 
in <h, in contrast to the interference I, which stems from 4)/ 

Remark. For the PPP, the EFIR does not depend on A, since Eq(/^^) oc 1/A. To see this, let 4>' = c<l> be a scaled 
version of <h. Then 

Ic = ^ /l 3 ;||x||““ = C““ ^ /l 3 ;||x||““ 
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and thus Ic^ = c^I ^■ Multiplying by the intensities, Ac^c ^ = XI ^ since X/Xc = c^- The same argument applies 
to all point processes for which changing the intensity by a factor c~^ is equivalent in distribution to scaling the 
process by c, i.e., for point processes <1>(A) where c<l>(l) = <h(c“^). This excludes hard-core processes with fixed 
hard-core distance but includes lattices and hard-core processes whose hard-core distance scales with 


Lemma 6 (EFIR for the PPP). For the PPP, with arbitrary fading, 

EFIRppp = (sinc(5)^/^. (36) 

Proof: The term in (l35l) can be calculated by taking the expectation of the following identity which 

follows from the definition of the gamma function r(x). 

Hence 

1 

= I i:}o,iJs)s-^+^ds. (37) 

From Slivnyak’s theorem ifTTl Thm. 8.10], EJ, = E for the PPP, so we can replace j (s) by the unconditioned 
Faplace transform Ci^{s), which is well known for the PPP and given by |[T6l 


= exp(—A7rE(/i^)r(l — 6)s^). 


From (ITTl) . we have 

1 

E{lJ) = 

r(o) Jq 

1 sine 5 

^ A7rE(/i^)r(l - (i)r(l + 5) ^ A7rE(/i^)' 

So X7r¥}g{If^)¥,{h^) = sine <5, and the result follows. ■ 

Remarkably, EFIRppp only depends on the path loss exponent. It can be closely approximated by EFIRppp « 
1-5. 


B. The tail of the SIR distribution 

Next we use the representation in (fT^ to analyze the tail asymptotics of the cedf Fsir of the SIR (or, equivalently, 
the success probability ps)- 

Theorem 4. For all simple stationary BS point processes <h, where the typical user is served by the nearest BS, 

Proof: From (|9l), we have Ps{6) = ¥,Fh{9R°‘I). Using the representation given in ([T^ . it follows from the 
Campbell-Mecke theorem that the success probability equals 




6l||x||" ^ hy\\y\\ 

y&^\{x} 


l(^>(6(o, ||x||)) = 0 ) 


= A /e’ 

FhioWxW^ ll^ll) einpty) 

J 

R2 

^ y6<I>x 
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where <I>a: — ^: y + ic} is a translated version of d>. Substituting x0^l^ 

!->■ X, 


Psi9) = xe 


-& 


I 




FhiWxr E ^y\ 


1 ( 6 ( 0 , ||x ||0 empty) 


dx 




= X9-^E: 


I E|,F (||x||"/oo) dx, 6^00 

Uij) [ hiWxDdx, e^oo, 


(38) 


where (a) follows since 9 0 and hence l{ 6 (o, ||x||0 empty} —)• 1. The equality in ( 6 ) follows by using 

the substitution x/^/“ —)• x. Changing into polar coordinates, the integral can be written as 


I (||x||") dx =7r6 / ^Fh{r)dr''= TrK{h'^), 

Jm Jo 

where (a) follows since h > 0 ifTTl . Since E(h) = 1 and 6 < 1 , it follows that E(h^) < 00 . 
For Rayleigh fading, from the definition of the success probability and Theorem HI 

-5 


(a) 


Ps(9) — ^isr(9) 


EFIR 


0 —>■ 00 . 


Hence the Laplace transform of ISR behaves as 0(9 for large 9. Hence using the Tauberian theorem in ifT^ 
page 445], we can infer that 

EFIR-^ 


E(ISR < x) ~ x'’ 

From Theorem IH the gain Goo immediately follows. 


r(i + 6)’ 


X ^ 0 . 


(39) 


Corollary 5 (Asymptotic gain at 0 —>• 00 ). For an arbitrary simple stationary point process 4> with EFIR given in 
Def HI the asymptotic gain at 9 ^ 00 relative to the PPP is 


Goo = 


EFIR 


/A7rEU7^)]^y/''^ 

\ sine 8 J 


EFIRppp 

Proof: From Theorem |4l we have that the constant c in ([ 8 ]) is given by c = EFIR^. cppp follows from Lemma 
[^as Cppp = EFIRppp = sine 6 . ■ 

The Laplace transform of the interference in (lT7l) for general point processes can be expressed as 


C: 


'o,iJs) = e:, (e 


-^Exe4.^»bll 


E[llCh{s\\x\\-^) = gl[Chis\ 




where is the probability generating functional with respect to the reduced Palm measure and Ch is the Laplace 
transform of the fading distribution. 

Corollary 6 (Rayleigh fading). With Rayleigh fading, the expected fading-to-interference ratio simplifies to 


where 


Ef\R=(\[ g-[A(x,-)]dx 

V Jr^ 

1 

A(x,y) = 


1/5 


1 + \\x\ 


a U, -Q 
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Proof: With Rayleigh fading, the power fading coefficients are exponential, i.e., Fh{x) = exp(—x). From (1381 ). 
we have 


p,{e)^xe-^ [ F}^F{\\xri)dx 


= xe 


-5 


<nT 




+ lkll"||y|| 


-dx, 


and the result follows from the definition of the reduced prohahility generating functional. 
For Rayleigh fading, the fact that 9^ps{9) sincJ as 6* —oo was derived in H] Thm. 2]. 


C. Tail of received signal strength 

While Theorem |4] shows that Ps{9) = &{9~^),9 oo, it is not clear, if the scaling is mainly contrihuted hy 
the received signal strength or the interference. Intuitively, since an infinite network is considered, the event of the 
interference being small is negligible and hence for large 9, the event S/I > 9 is mainly determined by the random 
variable S. This is in fact true as is shown in the next lemma. 

Lemma 7. For all stationary point processes and arbitrary fading, the tail of the ccdf of the desired signal strength 
S is 

P(5 > 0) ~ A7rE(/i^)0-^ 9^00. 

Proof: The cdf of the distance R to the nearest BS is Fr{x) ~ Attx^ for all stationary point processes ifTOl . 
Hence 

¥{S >9)= ¥{R < {h/9)^/^) 

~ X 7 rE [{ h / 9 )% 


So the tail of the received signal power S is of the same order Q{9 
affect the pre-constant. In the Poisson case with Rayleigh fading. 


and the interference and the fading only 


Ps(^) ~ A7rr(l + (5)0 0 —>■ oo. 

The same holds near 0 = 0. If for the fading cdf, FX(x) ~ ax"*, x ^ 0, 

F{S <9) = EFh{9R^) ~ a0”*E(i?”*“), 9 

For the PPP, 

FM -t- mrv /9.) 

0 ^ 0 . 


0 . 


P(5 < 0) ^ r(l + 'ma/2) 


(A7r)”*“/2 

So on both ends of the SIR distribution, the interference only ajfects the pre-constant. 

We now explore the tail of the distribution to the maximum SIR seen by the typical user for exponential h. 
Assume that the typical user connects to the BS that provides the instantaneously strongest SIR (as opposed to the 
strongest SIR on average as before). Also assume that 9 > 1. Let SIR(x) denote the SIR between the BS at x and 
the user at the origin. Then 

P(maxSIR(x) > 0) = E P(SIR(x) > 0) 

xe# 




= X 
= X 
= X9 


’^|,(SIR(x) > 0)dx 
1 


e- 


1 + 1 


X 


dx 


-5 


Go[Mx,-)]dx. 


From the above we observe that (for exponential fading), 

Ps(0) ~ P(niaxSIR(x) > 0), 0 —)• oo. 

xG# 

which shows that the tail with the maximum SIR connectivity coincides with the nearest neighbor connectivity. 
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Square lattice with Rayleigh fading fora=4 



Fig. 6. Scaled success probabitity pa(6)6^ for the square lattice point process with Rayleigh fading and a = 4. The asymptote (dashed 
line) is VEFIR ~ 1.19 and is tight as 0 —>■ oo. 


VI. Examples 


A. Lattices 


Let ui,U 2 be iid uniform random variables in [0,1]. The unit intensity (square) lattice point process $ is defined 
as For this lattice, with Rayleigh fading, the Laplace transform of the interference is bounded 

as l[T9l 


-sZ(2/5) < 

^ — ■^ 0 , 7 , 


(^)< 


1 

1 + sZ{2/6) ’ 


(40) 


where Z{x) = 4(^(x/2)/3(x/2) is the Epstein zeta function, Q{x) is the Riemann zeta function, and /3(x) is the 
Dirichlet beta function. Hence from (lT7l) 


Z{2/5)-^ < < 


TT csc(7r(i) 

T{5)Z{2/5Y' 


The upper bound equals {Z{2/5Yt{1 + <5) sinc<5) and it follows that for Rayleigh fading. 


Z{2/5) - Vsincii/ Z{2/5) 


(41) 


As a increases (5 —)• 0), the upper and lower bounds approach each other and thus both bounds get tight. 

The success probability multiplied by 6 ^, the EFIR asymptote and its bounds (|4T1) for a square lattice process 
are plotted in Figure |6] for a = 4. We observe that the lower bound, which is 1.29, is indeed a good approximation 
to the numerically obtained value EFIR « 1.40, and that for 0 > 15 dB, the ccdf is already quite close to the 
asymptote. 

For the square and triangular lattices. Fig. |7] shows the gain as a function of 6 and the asymptotic gains Gq 
and Goo for Rayleigh fading. Interestingly, the behavior of the gap is not monotone. It decreases first and then 
(re)increases to Goo- It appears that G{6) < max{Go,Goo}- If this holds in general, a shift by the maximum of 
the two asymptotic gains always results in an upper bound on the SIR ccdf. 

Fig. m shows the dependence of Gq and Goo on a. As pointed out in Subs. IIV-DI Gq is very insensitive to a. 
Goo appears to increase slightly and linearly with a in this range. 
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Horizontal gap for lattices, a=3 Horizontal gap for lattices, a=4 




Fig. 7. Gains G{9) for the square and triangular lattices and asymptotic gains Go and Goo (dashed) for Rayleigh fading and a = 3 and 
a = 4. 


Gg and for square and triangular lattices 



Fig. 8. Asymptotic gains Go and Goo (linear scale) for square and triangular lattices for Rayleigh fading as a function of a. 


B. Determinantal point processes 

Determinantal (fermion) point processes (DPPs) ll20l exhibit repulsion and thus can be used to model the fact 
that BSs have a minimum separation. The kernel of the DPP <1> is denoted by K{x, y) and—due to stationarity—is 
of the form K{x — y). Its determinants yield the product densities of the DPP, hence the name. The reduced Palm 
measure pertaining to a DPP with kernel is defined as 




1 


K{Xo,Xo) 


det 


K{x,y) K{x,Xo) 
K{xo,y) K{xo,Xo) 


(42) 
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Ginibre point process with Rayleigh fading, a=4 



Fig. 9. Scaled success probabitity ps{0)6^ for the GPP for Rayleigh fading with a 


4. The asymptote (dashed line) is at VEFIR « 0.89. 


whenever K{xo, Xq) > 0. Let K°{x, y) denote the kernel associated with the reduced Palm distrihution of the DPP 
process. The reduced prohahility generating functional for a DPP is given hy Il20ll 




n 


detf(l - (1 - /)iT°), 


(43) 


where detf is the Fredholm determinant and 1 is the identity operator. The next lemma characterizes the EFIR a 
general DPP with Rayleigh fading. 


Lemma 8. When the BSs are distributed as a stationary DPP, the EFIR with Rayleigh fading is 


EFIR = 


A f detf(l — (1 — A(x, •))iT°)da 
dR2 


1 /s 


(44) 


Proof: Follows from Corollary and (l43l ). ■ 

Ginibre point processes: Ginihre point processes (GPPs) are determinantal point processes with density \ = cj-K 
and kernel 

K{x,y) = 


Using the properties of GPPs lEH, it can he shown that 

OO n _ 

E!,(e-^^-) = n / 

7. 1 4o ^ 


Ig—cr 


k=l' 


^T{k) 


dr, 


from which Eq(/“'^) can he evaluated using (iTTl ). In Fig. |9j the scaled success prohahility 0^ps{9) and the asymptote 
EFIR'^ are plotted as a function of 6 for the GPP. We observe a close match even for modest values of 6 . Fig. [TO] 
shows the simulated values of the gains Go and Goo for the GPP as a function of the path loss exponent a. (^o ~ 1.5 
for all values of a, while Goo a j2. 


VII. Conclusions 

This paper established that the asymptotics of the SIR ccdf (or success probability) for arbitrary stationary cellular 
models are of the form 

p^{e) ~ 1 - co(9™, 0 0 ; p^{Q) ~ Coo^"^, 6 » oo 


for a fading cdf Fh{x) = 0(x™), x ^ 0. Both constants cq and Coo depend on the path loss exponent and the 
point process model, and cq also depends on the fading statistics. Depending on the point process fading may also 
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Gg and for GPP 



Fig. 10. Simulated gains Go and Goo for the GPP with Rayleigh fading as a function of a. 


affect Coo- Co is related to the mean interference-to-signal-ratio (MISR). For m = 1, cq = MISR, and for m > 1, 
Co depends on the generalized MISR. Cqo is related to the expected fading-to-interference ratio (EFIR) through 
Coo = EFIR"^. For the PPP, Coo = sine 5. The study of the MISR is enabled hy the relative distance process, which 
is a novel type of point process that fully captures the SIR statistics. A comparison of Go and Goo shows that a 
horizontal shift of the SIR distribution of the PPP by Go provides an excellent approximation of the entire SIR 
distribution of an arbitrary stationary point process. 

For all the point process models investigated so far (which were all repulsive and thus more regular than the PPP), 
the gains relative to the PPP are between 0 and about 4 dB, so the shifts are relatively modest. Higher gains can be 
achieved using advanced transmission techniques, including adaptive frequency reuse, BS cooperation, MIMO, or 
interference cancellation. As long as the diversity gain of the network architecture is known and the (generalized) 
MISR can be calculated (or simulated), the ASAPPP method can be applied to arbitrary cellular architectures. Such 
extensions will be considered in future work. A generalization to heterogeneous networks (HetNets) is proposed 
in |[22l . i^ . The method can be expected to be applicable whenever the MISR is finite. This excludes networks 
where interferers can be arbitrarily close to the receiver under consideration while the intended transmitter is further 
away, such as Poisson bipolar networks. 
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